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^ ■ Abstract 

Starch solutions, which are strongly non-Newtonian, show a surface instability, when 
subjected to a load. A droplet of the fluid is sandwiched between two glass plates and 
\ a weight varying from 1 to 5 kgs. is placed on the top plate. The area of contact 

' between the fluid and plate increases in an oscillatory manner, unlike Newtonian 

fluids in a similar situation. The periphery moreover, develops a viscous fingering 
like instability, which is not expected under compression. We attempt to model 
! the non-Newtonian nature of the fluid through a visco-elastic model incorporating 

generalized calculus. This is shown to exhibit a qualitatively similar oscillatory 
£q ' variation in the surface strain. 
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Forced spreading of a fluid under an impressed force is an interesting topic [U [2j 
[3] . The physics involved is challenging and the problem has important applications 
in technology as well. In real life many fluids are non-Newtonian, adding further 
complexity to the problem. 

In the present work we report studies on spreading of a starch solution between 
two glass plates, when the upper plate is loaded by a weight. We observe an in- 
teresting oscillation in the area of contact between the fluid and glass plate as a 
function of time. Earlier study of Newtonian fluids [3] did not show such behavior. 
We try to explain this phenomenon using fractional calculus, which is known to be 
an appropriate technique to study non- Newtonian, visco-elastic materials [IJ|5j|6]. 
Another remarkable observation is the appearance of a surface instability similar to 
viscous fingering. Viscous fingering under the condition of lifting, i.e. separating 
the plates is a well studied phenomenon [7J [SJ [§]. In this case the pressure is lower 
within the fluid, compared to the air pressure outside, satisfying the Saffman- Taylor 
condition for instability [10]. However, in the present case the fingering develops 
during compression. 

The fluid under study is prepared by adding m gms of arrowroot to 100 ml 
distilled water and boiling for 1 minute making the solution thicken. Here m = 



2.5 and 3.5 gms. A pinch of dye is added to enhance contrast. A drop of solution 
is placed on a glass plate and another glass plate is placed on top. The upper 
plate is loaded by a weight W, which varies from 1 to 5 kgs. The spreading drop is 
photographed from below using a video camera. The fractional change in the area of 
contact of fluid and glass is measured using Image-Pro Plus software and plotted as a 
function of time in Figured] A similar figure for ethylene glycol, which is Newtonian 
is shown for comparison in Figure |2j An oscillatory behavior is clearly visible for the 
starch solution. For the load of 1kg. the variation is nearly smooth, but oscillations 
are quite pronounced for 4 and 5 kgs. as shown more clearly in the inset. It is to 
be noted that this is not simply stick-slip behavior, where the strain would always 
increase, but in jumps. In this case, the strain actually decreases, before increasing 
again. To establish that the film does shrink and expand periodically, we show 
a superposition of two snapshots in figure© taken at an interval of 1 sec. Here 
the larger outer boundary corresponds to an instant of time earlier than the inner. 
So the ups and downs in figure© are genuine and not due to measurement error. 
Further, the corrugated appearance of the boundary demonstrates the instability 
developed. 

Rheological study of the two starch solutions was done at the Central Glass 
and Ceramic Research Institute, Kolkata by a Bohlin rheometer. Viscosity plotted 
against the rate of strain in Figure 01]) reveals the non-Newtonian nature of the 
fluids, which is more prominent for the higher concentration, as shown in the inset. 
This is in agreement with earlier results pTj . 

To analyze the stress-strain behavior here, we use the standard Maxwell model 
[12] of a spring and dashpot in series. The viscous term represented by the dashpot 
is further generalized by taking the gth derivative of the strain. Here q may be a 
fraction, for case q = 1 we have a Newtonian viscosity. The elastic term is assumed 
to be Hookean at present. The stress is then given by 



Here e is the strain, E the elastic modulus and f3 a parameter characterizing the 
effective viscosity of the non-Newtonian fluid. 

We assume a step function to represent the loading 



a(t) = a for t > 0, and a(t) = for t < 



. The initial condition for strain is e(t) = for t < 0. 
The Laplace transform of equation© gives 



w E [s st + E/fi 
The inverse Laplace transform of equation© gives 
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where, ML(-kt) is the one parameter Mittag-LefHer function defined by 



ML^-kt) = e~ kt 

oo 

ML q (z) = 



We now plot the strain as a function of time using equationQ. For q = 1, the 
strain increases smoothly and saturates to the value = a/E. For q < 1, a similar 
behavior is observed, but the saturation value is less. For q > 1 however, we see an 
initial increase in strain overshooting followed by oscillations before saturating 
to eoo. The oscillations are more pronounced as q increases. These results are shown 
in figure ()5]). 

The equation [1] is the generalized representation of the stress-strain relation. Let 
us rewrite the same as: 

"IF + Be(t) = r (t) (4) 

Where B = E//3. Here (3 is the generalized viscosity coefficient with units cor- 
responding to the non-integer order. When the order q — 1, then normal coefficient 
of viscosity is recovered. The unit of B for order q — 1 is per seconds i.e. s _1 , but 
for any other order q ^ 1; the unit modifies to s~ q . 

Mathematically one has to see the Green's function for general relaxation in 
equatiodT], so we write the homogeneous equation with RHS equal to zero. To 
that, we give delta function stress excitation. The strain built up for any relaxation 
process may be treated as convolution integral of a strain variable with integral 
kernel K q (t), as 0(6]. 

j t <t) = - jT K q (t - r)e(t)dt (5) 

Well if the memory kernel is K(t) = B 5(t), we have the above system [T] without 
memory [51 16] and the Green's function will be 



e(t) = e e 



-B t 



that is the impulse response quickly decays to zero. Here e is initial strain of the 
system at t = 0. 

This can be derived as follows: 

K(t) = B 8(t) (6) 

d 



f 6(t - r)e(t)dt = -B e(t) (7) 
Jo 

e(t) = e e- B °* (8) 



The homogeneous strain relaxation equation for no-memory case is a first order 
Ordinary Differential Equation i.e. 



j t e{t) + B e(t) = (9) 

If the memory kernel is a constant say K2(t) = B 2 , then we will have oscillatory 
Green's function, which never decays to zero [5], [6] . 

K 2 (t) = B 2 (10) 

^e(t) = -B 2 e(t) (11) 

e(t) = e cos(^ 2 t) (12) 
The generalized memory integral is [51 E] . 

K{t) = B q t q - 2 ; < q < 2 (13) 



d . . 1 

— e (t) = 

dt w ri 



(14) 



r q = [B q T(q-l)]- L (15) 

Its corresponding generalized differential equation, obtained from above deriva- 
tion, is the system with memory with the memory index coming as fractional order 
of the Fractional Differential Equation with, < q < 2. 

d q e(t) t~ q „ , , 

— - e — = -r~ q e(t) 16 

difl T(l-g) w v ; 

In our experiment, to describe the oscillatory response to a step input we say that 
the order is between 1 < q < 2 and thus the system has a long lingering memory. 
In equation f|T6|) above, if the initial stress be eo , assuming Heaviside's step function 
as the stress input, it modifies to equation([T]), for r~ q — > B. 

The order of the equation([T|) corresponds to a system with memory. The non- 
Newtonian fluids without oscillatory behavior will have < q < 1 , which is fractional 
order, and the step-response will have monotonically increasing strain response, 
given by one argument Mittag-Leffler function. Its impulse response will be having 
long tailed decay. That is the response will have long-range temporal correlation. 
The Newtonian fluid will have integer order in equation([T|) with q — 1, representing 
the system without memory, and the step-response will have monotonically increas- 
ing strain as e(t) ~ 1 — exp(Bot); where its impulse response will decay quickly as 
e(t) ~ exp(Bot). This visco-elastic system with Newtonian viscous bahavior can be 
modeled with a discrete ideal spring and a ideal dashpot. Whereas the more com- 
plicated case with non-Newtonian viscosity requires a different representation like a 
fractal chain of the ideal spring and ideal dashpot combination. We have observed 



oscillatory strain and thus infer the fractional order q of our system to lie between 
1 and 2. 

Considering that B = ^ is a parameter representing the relative strengths of the 
elastic and viscous terms in equation([T|), we may see how the system responds to 
changes in B. We find that variation of B changes the time period of the oscillations 
in strain, without affecting the amplitude as shown in figure(j6]). The time period is 
smaller when elasticity dominates. The amplitude depends on q and of course the 
magnitude of the strain changes proportionately to the load a, also affecting the 
amplitude of oscillation. 

If we assume that the fractional change in area is equivalent to the strain in 
the system, the results for the starch solutions are reproduced qualitatively by the 
visco-elastic model with non-Newtonian rheology. Here the qth order derivative 
takes care of the non-linearity in the complex fluid. In the earlier paper [3] the 
Newtonian fluids were assumed incompressible and the change in film thickness was 
calculated from the constant volume of the fluid. Here the fluid may have a finite 
compressibility, so we refer only to the area which is measured directly. 

To conclude, this work demonstrates the interesting phenomenon of oscillatory 
spreading of starch solutions on glass and illustrates further how the approach of 
generalized calculus may be used to analyze it. The qualitative agreement of the 
variation in strain with time from the experiments and the theory is encouraging 
and further exploration along these lines promises to be rewarding. 
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Figure 1: Experimental variation in strain with time for different values of load for 
2.5 concentration arrowroot solution on glass. The inset shows the oscillations at 
higher magnification. Symbols - diamond (black), plus (red), triangle (green), circle 
(blue) and square (violet) represent respectively the effects of loads 1,2,3,4 and 5 
kgs. in addition to the glass plate which weighs around 0.6 kgs. 




Figure 2: Variation in strain with time for different values of load for ethylene glycol 
on glass. 




Figure 3: A snapshot of the film (inner light colored blob) superposed on the pho- 
tograph of the film photographed 1 second earlier (darker outline visible along the 
periphery) shows the shrinking of the film. 




Figure 4: Variation of apparent viscosity with strain rate for the arrowroot solu- 
tions shows the non-Newtonian behavior, which is stronger for higher concentration 
(inset). 




Figure 5: Variation in strain with time for different values of q. 




Figure 6: Variation in strain with time for different values the parameter B for q 



